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A systematic study of both classical and quantum geometric frustrated Ising models with a com- 
peting ordering mechanism is reported in this paper. The ordering comes in the classical case from 
a coupling of 2D layers and in the quantum model from the quantum dynamics induced by a trans- 
verse field. By mapping the Ising models on a triangular lattice to elastic lattices of non-crossing 
strings, we derive an exact relation between the spin variables and the displacement field of the 
strings. Using this map both for the classical (2-|-l)D stacked model and the quantum frustrated 
2D system, we obtain a microscopic derivation of an effective Hamiltonian which was proposed 
before on phenomenological grounds within a Landau-Ginzburg- Wilson approach. In contrast to 
the latter approach, our derivation provides the coupling constants and hence the entire transverse 
field-versus-temperature phase diagram can be deduced, including the universality classes of both 
the quantum and the finite-temperature transitions. The structure of the ordered phase is obtained 
from a detailed entropy argument. We compare our predictions to recent simulations of the quan- 
tum system and find good agreement. We also analyze the connections to a dimer model on the 
hexagonal lattice and its height profile representation, providing a simple derivation of the contin- 
uum free energy and a physical explanation for the universality of the stiffness of the height profile 
for anisotropic couplings. 
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I. INTRODUCTION 

Frustrated two-dimensional magnetic systems exhibit 
a rich variety of phases and critical points. The conse- 
quences of quantum or thermal fluctuations about the 
highly degenerate ground states of frustrated magnets 
are important to understand in the ongoing quest to find 
novel exotic quantum statesi. The past decades have 
seen a resurgence of significant interest in a systematic 
study of geometrically frustrated antiferromagnets^. Ge- 
ometric frustration arises in materials containing anti- 
ferromagnetically coupled moments which reside on ge- 
ometrical units, that inhibit the formation of a coUinear 
magnetically ordered Neel state, and induces a macro- 
scopically large degeneracy of the classical ground state. 

For magnets with a discrete Ising symmetry, the com- 
plexity of the ground-state manifold can endow the sys- 
tem with a continuous symmetry. Such symmetry is of 
particular importance to 2D quantum magnets at low 
temperatures since then the Mermin- Wagner theorem ap- 
plies, precluding an ordered phased. Another possible but 
contrary scenario is "order-from-disorder"— where zero- 
point fluctuations select a small class of states from the 
ground-state manifold since those states are particularly 
susceptible to fluctuations. This fundamental mechanism 
can produce an ordered symmetry-broken stated. Hence 
one expects weak competing fluctuations about the clas- 
sical ground states to be able to generate new strongly 
correlated states and (quantum) phase transitions of un- 
expected universality classes. 

From a theoretical perspective, it is natural to try to 
understand the role of frustration from extremely sim- 
ple interactions and dynamics. The possibly simplest 
realization of classical frustration is found in the antifer- 
romagnetic Ising model on a triangular lattice (TIAF). 



Each elementary triangle is frustrated, and the TIAF is 
disordered even at zero temperature with a flnite entropy 
density and algebraic decaying spin correlations™. The ef- 
fect of quantum fluctuations about the highly degenerate 
ground states can be studied in its simplest form by in- 
troducing quantum spin dynamics from a magnetic field 
which is transverse to the spin coupling. For this trans- 
verse field TIAF model, and its companions on other 2D 
lattices, Moessner and Sondhi have argued the existence 
of both ordered and spin liquid phases^. In the limit of 
a small field, the quantum ground state is constructed 
as a linear superposition of classical ground states which 
maximize the number of spins which can be flipped to 
gain transversal fleld energy at no cost in exchange en- 
ergy. This yields a strong suppression of conflgurations 
and, since the TIAF is already critical at zero field, or- 
der emerges. As a result, there will be a discontinuity in 
entropy and correlations in the ground state at vanishing 
field and zero temperature. At a large transverse field, 
the TIAF has a gaped paramagnetic ground state. It is 
separated from the ordered state by a quantum critical 
point. The additional effect of thermal fluctuations has 
been studied quantitatively so far only in simulations^. 

Experimental realizations of these frustrated Ising sys- 
tems can be found either directly in magnets with strong 
anisotropy, e.g., in LiHoF^-, or indirectly (via the equiv- 
alent (2-|-l)D classical model) in stacked triangular lat- 
tice antiferromagnetsS with strong couplings along the 
stacking direction as studied in recent experiments on 
CsCoBra"'^". However, transverse field Ising models can 
also provide insight into some more complicated systems 
in certain limits. They may describe the singlet sector be- 
low the spin gap of frustrated antiferromagnetic quantum 
Heisenberg models, e.g., on the Kagome lattice, since the 
latter model can then be formulated as a Z2 gauge theory 
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which in turn is related by duahty to the transverse field 
Ising system on the dual lattice of the original Heisenberg 
model^^. An ordered Ising phase describes then in the 
Heisenberg problem a paramagnet with spin waves form- 
ing the gaped excitations. A related approach to study 
Heisenberg antiferromagnets is via quantum dimer mod- 
els (QDM) on the same lattice^-. Depending on the lat- 
tice symmetry, the QDM exhibits a disordered state (on 
the triangular lattice)i2*i^, corresponding to a resonating 
valence bond or spin-liquid phase, or it is found always 
in ordered valence-bond solid phases (on the hexagonal 
lattice)iii^. Interestingly, the spin liquid phase renders 
the triangular QDM a promising candidate for quantum 
computingi^. For Ising spins, there is also a direct cor- 
respondence between QDM's with just kinetic terms and 
fully frustrated Ising models in a small transverse field 
on the dual lattice^. In particular, the transverse field 
TIAF maps to the hexagonal QDM. 

For comparison to the new developments presented 
in this work, we briefly review the salient achievements 
for the classical stacked and quantum TIAF obtained in 
earlier studies. Symmetry arguments have been used 
to guess a Landau-Ginzburg- Wilson (LGW) theory for 
the (2-1- 1)D classical stacked system^^. After transla- 
tion to the quantum TIAF, this approach suggests a 
quantum critical point of 3D XY universality at inter- 
mediate transverse field strength and an extended criti- 
cal phase at finite temperatures^. However, the Suzuki- 
Trotter mapping relating the 2D quantum and (2-1- 1)D 
classical systems involves a scaling limit with a diverg- 
ing anisotropy of the classical exchange coupling which 
may invalidate a LGW approach^^. Moreover, the latter 
approach has been put somewhat into question, mainly 
since it has been argued that it fails to describe the (clas- 
sical) system at low temperatures since it neglects the 
restriction to classical spin values ±1 and thus geomet- 
rical frustrationi^. Very recent Monte Carlo simulation 
support the LGW based conjecture for the phase dia- 
gram but the actual computations were performed for 
the (2-1- 1)D classical problem^. Recently, the weak field 
behavior has been also studied in terms of a quantum 
kink crystali^. 

In a recent letter— the present authors have presented 
a string description of the transverse field TIAF in or- 
der to obtain a quantitative prediction for the phase di- 
agram. In the present paper, we will give a full account 
of the relation between spins and non-crossing strings, 
including a more complete presentation of the relation to 
dimer models and its surface height representation. We 
exploit and extend the quantum dimer analogy in order 
to map the TIAF, at arbitrary transverse field strength, 
to quantum strings which result from the superposition 
of dimer configurations plus defects and a fixed reference 
dimer state. Using the Suzuki- Trotter theoremii we ob- 
tain a stack of coupled 2D layers of classical strings which 
in the Suzuki- Trotter limit of infinitely many layers can 
be described by the Villain model2i. We show that the 
LGW actionpi^, which was employed in the vicinity of 



the quantum critical poin'b^, can be derived microscopi- 
cally from the quantum string action if the phase of the 
complex LGW order parameter is identified with the dis- 
placement field of the string lattice. However, our ap- 
proach explicitly takes into account frustration which is 
encoded in the topological constraint on the phase field 
resulting from the non-crossing property of the strings. 
This constraint restricts the phase field configurations of 
the order parameter and thus distinguishes our theory 
from the original LGW approach. Our approach con- 
firms explicitly the 3D XY universality of the quantum 
critical point, and allows to predict the phase diagram 
at arbitrary transverse field strength and finite tempera- 
ture. Using an entropy argument for the strings, we can 
determine the nature of the ordered phase. 

In this paper we present first a thorough discussion of 
the relations between Ising spins on the triangular lattice, 
dimers on the hexagonal lattice and its height represen- 
tation and elastic string lattices. After a introduction of 
the spin models in Sectional we explore the relations to 
dimers and strings in Sect ion imi for the classical 2D Ising 
model on the triangular lattice in order to set up the for- 
malism for the following sections. In Section [l VI we map 
the Ising model on a stacked triangular lattice to (2-|-l)D 
string lattice which is described by a 3D XY model with 
a 6-fold symmetry breaking term. Using the results of 
the latter section, we perform a Suzuki- Trotter mapping 
of the 2D quantum Ising system to the classical stacked 
model which enables us to predict a quantitative phase 
diagram for the quantum frustrated model in Section FVl 
Finally, we conclude with a summary and a discussion of 
potential extensions of our work in Section [VII 



II. MODELS 

We are interested in two models which are both based 
on the Ising antiferromagnet on a triangular lattice 
(TIAF). First, we study a classical three dimensional 
Ising system which consists of a stack of TIAF's which 
are coupled ferromagnetically. The Hamiltonian reads 

^^3D = •-'^ll ^ CTikCTjk — J±''^(^ik(^ik+1, (1) 

{ij),k i,k 

with aik = ±1 and where Jy, J± > and (ij) in- 
dicates summation over nearest-neighbor pairs in each 
TIAF plane. This system is fully frustrated in each 
layer, but has no competing interaction along the stack- 
ing direction. This model has been studied originally 
within a Landau-Ginzburg- Wilson (LGW) approach by 
Blankschtein et al}^. In fact, much of the physics of the 
stacked triangular lattice is similar to that of the two- 
dimensional triangular lattice. However, the presence of 
the third dimension has the tendency to stabilizes or- 
dered phases. Indeed, for all well-characterized mate- 
rials which have triangular magnetic lattices, the long- 
range order at low temperatures is three-dimensional in 
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nature. The stacked system is of direct experimental rel- 
evance since it describes the low-temperature physics of 
the Ising-like compounds CsCoCla, CsCoBra and related 
materials where transition metal atoms form chains along 
the stacking direction which are coupled through three 
equivalent halogen atoms. For a review on these com- 
pounds see Ref. '5. Almost the entire existing body of 
theory on this stacked system is based on the LGW ap- 
proach, mean-field theory and numerical simulations with 
sometimes conflicting results for the low-temperature or- 
dered phase. We will briefly review previous develop- 
ments for this model in Section Hvl 

As shall be demonstrated below, the previously in- 
troduced model is also useful in describing the two- 
dimensional quantum antiferromagnet which results from 
applying a transverse magnetic field to the TIAF. The 
latter system has the Hamiltonian 

i7 = afa|-|-r^<, (2) 

tj^, are Pauli operators, and F is the transverse field. 
This model is of particular interest since it combines a 
standard realization of geometric frustration with simple 
quantum dynamics. The transverse field induces tun- 
nelling between the exponentially large number of clas- 
sical ground states at zero temperature. One can argue 
that the quantum fluctuations select a smaller suscepti- 
ble class of the ground states. This would lead then to 
a reduction of ground state entropy, and an "order from 
disorder'"^ phenomenon is expected. 

Both the classical stacked magnet and the 2D quantum 
TIAF are related to the quantum dimer model (QDM) 
on the hexagonal lattice. 

III. STRING PICTURE OF CLASSICAL 
FRUSTRATED 2D MAGNETS 

A. From spins to dimers 

We start with the description of the antiferromagnetic 
Ising model on the triangular lattice in order to explain 
its relation to dimers and strings. The Hamiltonian is 

Htiaf = ^ Jij<Ti(Jj (3) 

{i.j) 

where the couplings Jij are equal to one of the three pos- 
itive constants Ji, J2 or J3 depending on the orientation 
of the bond (ij) relative to the three main lattice direc- 
tions. It is well known^ that the ground state of this sys- 
tem has exponentially large degeneracy in the isotropic 
case Ji = J2 — J3, leading to a finite entropy density 
and the absence of order even at zero temperature. For 
all ground states of the general model of Eq. Q each tri- 
angle has exactly one frustrated bond, and hence there is 
a one-to-one correspondence (modulo a global spin flip) 
of all spin ground state configurations to complete dimer 



FIG. 1: (color online) Mapping of a spin configuration to a 
dimer covering on the dual hexagonal lattice. 

coverings of the dual hexagonal lattice. This can be seen 
from Fig.^ If one places a dimer across every frustrated 
bond, one obtains the corresponding dimer covering of 
the dual hexagonal lattice, where each lattice site in the 
hexagonal lattice is touched by one and only one dimer^Si. 
If the two smallest couplings are equal, Ji — < J2, 
then dimers can occupy bonds in directions 1 and 3 only 
and there is still a huge ground state degeneracy. But the 
entropy at T = scales now ~ ^/N with the number of 
lattice sites, TV, due to the constraint that dimers cannot 
touch each other. Thus the entropy density vanishes and 
the system is ordered at T = but disordered at any 
finite temperature. (At finite T topological defects have 
to be considered, see below.) If one of the three coupHng 
constants is the smallest, then there is only one dimer 
ground state configuration with all dimers being perpen- 
dicular to the direction with smallest coupling, and order 
survives at finite temperature. 

The analogy to dimers can be expressed by the fact 
that the Ising partition function in the limit of vanishing 
temperature is proportional to that of the dimers which 
reads 

Zn = J2^i'^T, (4) 

{D} 

where the sum runs over all complete dimer coverings of 
the hexagonal lattice, and ni and 712 is the number of 
dimers being perpendicular to direction 1 and 2, respec- 
tively. Here we assumed that J3 < Ji, J2 so that the 
weights for the two types of non- vertical bonds are given 
by 

and all vertical bonds have a weight of unity. Of course, 
for the TIAF the isotropic case with Ji = J2 = J3 is of 
most interest but for the dimer model itself anisotropic 
weights provide some interesting insight as will become 
clear below. We note that even at finite temperature the 
analogy to dimers will turn out to be useful but then 
defects (triangles with all three bonds crossed by dimers) 
have to be included. 
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FIG. 2: (color online) Mapping a dimer covering (a) to a string 
configuration (c) via the reference state (b). The numbers de- 
note the height profiles associated with the dimer coverings. 
The displacement of a string from its straight reference posi- 
tion is determined by the height profile at the two plaquettes 
which are joined by the displaced dimer in the original dimer 
covering (a). 



B. From dimers to strings 

In this subsection we make use of another mapping 
which relates dimer configurations to fluctuating strings 
which are directed and non-crossingS^. This mapping 
appUes to general weights for the dimers but of primary 
interest is the TIAF at T = where the dimer weights 
are either unity or zero, and there is a one-to-one cor- 
respondence (modulo a global spin flip) of Ising ground 
states to string configurations. The string representation 
results from the superposition of a given dimer covering 
and a fixed reference dimer covering where all vertical 
bonds of the hexagonal lattice are occupied by dimers, 
see Fig.|2Ib). The superposition is an "exclusive or" op- 
eration, i.e., only if a given bond is covered by a dimer 
either in the original covering or in the reference cover- 
ing, it will be covered in the superposition. The resulting 
superposition is no longer a dimer covering (since dimers 
touch each other) but an array of strings which are di- 
rected (along the reference direction) and non-crossing 
(due to the fact that each site is touched by exactly one 
dimer in the original dimer covering)24. The fluctuations 
of the lines result from the non-zero entropy density of 
the spin system. Or in other words, we map the zero 
temperature TIAF to a line lattice at a finite (virtual) 
temperature Tj. 

Before we parameterize the string configurations, it is 
useful to introduce a height profile h which can be as- 
sociated with a dimer covering. If the lattice covered 
by dimers is bipartite, then there exists a well-defined 
(single- valued) height profile which is defined on the sites 
of the original (triangular) lattice. Starting at an arbi- 
trary site with some integer number, one follows a trian- 
gle pointing down clockwise and changes the height by 
+2 (-1) if a (no) dimer is crossed. Repeating the latter 
process for all triangles pointing downward, one obtains 
a consistent height on all sites, see Fig. |21 (a). If one 
subtracts the height profile of the fixed reference dimer 
covering from a given height profile, one observes that the 
previously introduced strings separate domains of equal 



FIG. 3: Different kind of steps of the random walk of a string. 

height which is an integer multiple of 3, see Fig. |2fc). It 
will turn out that the height profile measures the rough- 
ness, i.e., the displacement of the strings from a perfectly 
straight configuration. 

The array of strings is characterized only by its density 
p ~ l/£ and the elastic line tension g of the strings. The 
interaction between the strings consists only of the non- 
crossing constraint and thus introduces no additional en- 
ergy scale. The configurations of a single directed string 
are characterized by its position function Xi{z) since over- 
hangs are forbidden. The reduced energy of all strings is 
then 

First, we determine the line tension g from a simple ran- 
dom walk argument. If x denotes the mean string po- 
sition, then the total mean squared displacement of a 
string of length L = y/3aM along the y-direction (after 
M steps), with a the triangular lattice constant, must be 

M 

Y,{{x,-xf)=m/g (7) 

i=l 

in order to be consistent with the result of the contin- 
uum model of Eq.@. Here x^ = —a, 0, +a are the 
three possible transversal steps. The different kind of 
steps correspond to the paths OL, OC and OR shown 
in Fig. 131 The corresponding probabilities for the steps 
can be expressed in terms of the weights of the occupied 
non-vertical bonds as 

Z2 2z\Z2 z\ 

" {z^^Z^r " (^1+^2)2' " (^1+^2)2- 

(8) 

The mean position after one step is 

X ^ {pB^ - Pl) a ^ ^^—^a. (9) 

Zl + Z2 

Thus, after M steps, the variance of the transverse wan- 
dering around the mean position is 

M 

Y^iix. - xf) ^ L a, (10) 
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where rj = zi/z2, yielding the Hne tension 
g V3(2 + 7? + 7r') 



Ti 



2a 



(11) 



The mean density of strings depends also on the dimer 
weights of the non- vertical bonds. It can be computed 
exactly with the resul t - ?'^^ 



P - 



{Zl + Z2f 



1 



4ziZ2 



1/2 



(12) 



for Z\ -\- zi > 1. For zi + Z2 < 1 the den- 
sity vanishesS^, corresponding to a Kasteleyn or 
commensurate-incommensurate transition. For the 
isotropic TIAF, one has zi = Z2 = 1 and p = 2/ (3a). 

Before we introduce the relation between the string 
displacement and the height profile, and a correspond- 
ing continuum description, it is instructive to focus on 
the free energy of the string array. This will allow for 
an independent check of our random walk description. 
Moreover, from the free energy one can determine the 
elastic constants on large length scales which are needed 
for the proper continuum model. 

Actually, the density of strings is determined by the 
densities pi, p2 of non- vertical dimers with weights zi, 
Z2, respectively. This can be seen from the fact that 
along a cut perpendicular to the string direction each 
non-vertical dimer corresponds to one string. It follows 
from Eq. (gj that 



Pi 



Zj d log Zp 

N dz. 



(13) 



for i = 1, 2, where the number of triangular lattice sites, 
N, equals the total number of dimers. The string density 
then reads 



P = 



Pi + P2 



(14) 



Now we change variables from zi, Z2 to 77 = z\l z^ and 
p = pa, see Ea. lfT^ . Since zid'q/dzi-\-Z2drj/dz2 = 0, and 
using Ea. lfT^ . we have 



9 ^ d I + rf + 2t] cos{tt p) d 
dzi dz2 TTi] sin(7rp) dp 



(15) 



This allows us to derive a differential equation for the 
free energy density of the dimers. 
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(16) 



as a function of the parameters ?/, p of the string system. 
fo fulfills the equation 



dfo 27r 77psin(7rap) 



dp a/3 l + if + 27] cos{Trap) ' 



(17) 



Integrating the latter expression yields the exact result 
for the lattice string model. In the continuum limit a — > 
0, the result simplifies to 
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D 



3 V3(2 + r; + 77-i) 



(18) 



There is an alternative and simple approach to obtain 
the free energy in the continuum limit. The strings can 
be regarded as the world lines of free fermions in one 
dimensions^. The non-crossing condition is then imple- 
mented automatically by the Pauli principle. In the limit 
of infinitely long strings, the ground state energy den- 
sity of the fermions equals the free energy density of the 
strings if the fermion mass is mapped to the line tension 
g and h to the temperature T;. From the expression for 
the ground state energy of one-dimensional free fermions 
follows immediately the reduced free energy density of 
the strings. 
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(19) 



where /i is the reduced free energy (per length) of a single 
string, and the last term is generated by the reduction 
of entropy due to the non-crossing condition. A direct 
comparison of the free energies of Eas. l|18|l and H19|) is 
not possible since the number of strings varies with the 
dimer coverings although the number of dimers is fixed. 
Only the mean density of strings is fixed. Hence, one has 
to consider the grand potential density j for the string 
system which follows easily from the free energy density. 



3 = Ji- PP= — z P ' 

3 9 



(20) 



where fi — dfi/dp is the chemical potential (per string 
length). Comparing Ea. (|18|l and Ea. l|2U|) . we see that 
fo = j exactly for the expression of the line tension 
which we obtained above from an independent random 
walk argument, see Ea. (|lll) . Thus in the continuum limit 
the ground states of the TIAF and the dimer model can 
be described as free fermions with their mass determined 
by our simple random walk argument on the lattice. 

Now we turn to the structural properties of the string 
lattice. The most ordered state consists of straight 
("flat") strings which do not wander transversally. De- 
pending on the string density there are different com- 
mensurate and incommensurate states possible. For the 
discussion of the flat states we focus on the state with 
a maximal density of p = 2/ (3a) which corresponds to 
the most interesting case of isotropic couplings for the 
TIAF. In this case there are two non-equivalent classes 
of flat strings which are not related by shifts by a nearest- 
neighbor vector of the triangular lattice, see Fig.^ These 
classes translate directly to different ordered spin states 
which are characterized by their sublattice magnetiza- 
tions. One state, (-1 — 0), is obtained by orienting two 
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(+-0) (++-) 

FIG. 4: Two difFerent flat states, with the corresponding sub- 
lattice magnetizations of the isotropic TIAF indicated below. 

sublattices uniformly but opposite while the third sub- 
lattice has the same number of -I- and — spins, leading to 
zero magnetization. This, on average, locks the straight 
strings on the sublattice with zero magnetization. Since 
each of the three sublattices can be chosen as the one 
with zero magnetization, equivalent flat string states are 
related by a shift of a/2. The other state, — ), has up 
spins on two sublattices and down spins on the other sub- 
lattice which places the non-vertical dimers on straight 
strings that are locked symmetrically between the two 
sublattices with equal magnetization. 

The displacement u from flat strings is determined by 
the original height profile h of the corresponding dimer 
covering. Since only non-vertical dimers are conserved 
under the subtraction of the reference covering, we will 
defined u for each non-vertical bond of the hexagonal 
lattice. If a bond is occupied by a string segment, then 
u measures the normal distance between the correspond- 
ing flat string and the center of the segment, cf. Fig.|3c). 
Since a shift of the height hhy 3 corresponds to a string 
translation by the mean string distance the displace- 
ment on each non- vertical bond is given by u = £h/3 + uo 
where h is the coincident (original) height of the two 
hexagonal plaquettes which are joined by the non- vertical 
bond, see Fig.[2Ia). uo is a global offset of the flat string 
lattice. For uq — na/2, n G Z, the class — ) is chosen 
while Uq — a/A + na/2 selects (-1 — 0). This definition of 
u is equivalent to the height profile introduced by Zeng 
and Henleji21 by averaging over the three sites of every 
triangle in order to obtain a coarse-grained height on the 
center of the triangles. 

After a coarse-graining over length scales large com- 
pared to the lattice constant, one obtains a continuous 
field u{r) which allows to write the effective free energy 
of long- wavelength fiuctuations of the string lattice in the 
form of a continuum elastic energy, 

ld'r{'f{d^u)' + 'f{dyur + V,(u)} (21) 

with compression cn and tilt C44 modulus. The tilt mod- 
ulus is fixed by the line tension g as C44 = {g/Tijp. 



The compression modulus on asymptotically large length 
scales is affected by the entropic repulsion between 
strings and has to be determined in a macroscopic way. 
The compressibility of the lattice is given by the second 
derivative of the free energy density with respect to the 
mean distance i — 1/ p between the stringsS^, which leads 
to the result 

c,,^i^[m)]^^'''^p\ (22) 

where we have used Eq. p9|l . Vl (u) is a periodic potential 
which reflects the discreteness of the lattice and favors 
an ordered phase with flat strings. Since equivalent flat 
states are related by shifts of all straight strings by a/2 
and since u = th/2> -\- Uq, the locking potential must have 
the formes 

Vl — —vcos — (u — uo) = — ucos \ — — h] (23) 
V a J V 3 a / 

with V > 0. If Y^ciiC44 > 27r/a^ the string lattice is 
sufficiently stiff and the potential Vl is relevant (under 
renormalization) so that the strings lock into one of the 
fiat states. Interestingly, y^cnciJ = n/i'^ is independent 
of the line tension, and the locking potential is relevant 
for i < a/\/2. However, the minimal string separation is 
£ = 3a/2 (corresponding to isotropic dimer weights), and 
hence Vl is irrelevant at all possible densities. 

However, if there are additional interactions present 
which increase the stiffness of the string lattice, the lock- 
in potential might become relevant. In fact, we will see 
below that a coupling of many TIAF can render the pe- 
riodic pinning relevant. Therefore, it is important to de- 
termine the consequences of the lock-in potential for the 
spin configurations, which depend on the sign of that po- 
tential. To do so, we recall that the offset uq selects the 
lattice positions of the strings in the flat states. It has 
to be determined from additional information which is in 
the present case provided by spin configurations which 
correspond to flat strings and small fluctuations around 
the locked-in states. Thus, we focus now on the isotropic 
TIAF so that the possible classes of fiat states are those 
shown in Fig. 0] If one switches between the two classes, 
Uo shifts by a/4 and thus the sign of Vl changes. The 
sign can be determined as follows. 

The systems selects the class with the largest en- 
tropy, i.e., with the maximal number of configurations 
which yield at arbitrary large r still a finite displacement 
([u(r) — u(0)]^)flat ~ const., where the displacement is 
measured relative to the respective flat state. For large 
systems, the entropy of such macroscopically fiat states 
can be estimated from the number of string configura- 
tions on the hexagonal lattice for which each string fluc- 
tuates in a tube of width i about its straight reference 
position, i.e., |u| < £/2 = 3a/4 for each non-vertical 
dimer—. This puts a strong constraint on the spin states 
since a flip of a single fiippable spin (a spin with 3 up 
and 3 down spins as neighbors) shifts a string by the 
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hexagonal plaquette on which the spin sits. In the class 
(H — 0) the spins on one sublattice (dark grey in Fig. ^ 
can assume any of the 2^/^ possible states where N is the 
number of triangular lattice sites. The maximal displace- 
ment of u = 3a/4 is obtained if spins on the other two 
sublattices are also flipped with respect to the perfectly 
ordered state (+ H — ) which, however, is possible only 
for 1/2^ X 2N/3 sites since the 3 neighbors on the fully 
flippable sublattice must have the same orientation as 
the spin to be flipped in order to have directed and non- 
crossing strings. This yields 2^/^+^/12 = cxp(0.2888 TV) 
configurations. For the fluctuations in class (-1- H — ) all 
spins on one sublattice (white in Fig. 2)) are frozen. The 
sites of the other two sublattices can be divided into iV/12 
"rings" , each of 6 plaquettes centered about a frozen site, 
and extra plaquettes between the rings. Due to the 
constraints on the strings, each ring permits only 18 spin 
states. For a given spin state on all rings, the extra pla- 
quettes can only flip if its 3 neighboring spins located on 
rings point up, which occurs with probability (13/18)'^. 
Thus there are IS^/^^ x 2(i3/i8)'w/6 = exp(0.2844 TV) 
conflgurations. We conclude that entropy favors flat 
states of the class (H — 0) and one has to set uq = —a/4 
as shown in Fig. |21 (c), i.e., Vl{u) = ucos(47ru/a) with 
positive amplitude v. 



C. Spin-spin correlations 

In this section, we are going to apply the above frame- 
work to the spin correlations of the TIAF model at zero 
temperature. We concentrate on isotropic couplings since 
otherwise the system is ordered. We have seen that 
the ground state manifold of the isotropic TIAF can be 
regarded as the configuration space of a string lattice. 
Based on the irrelevance of lattice pinning effects, one can 
expect from the effective elastic description of Ea. H21l) an 
algebraic decay of correlations which is characteristic for 
such kind of two-dimensional systems. In fact, this cor- 
responds to the exact result found by Stephenson^^. 

The displacement correlations of the string lattice can 
be easily computed for large r from Ea. (|21|l with the 
result 



i^logf) 

\/CllC44 \a/ 



(24) 



So far, we have constructed the string displacement field 
and the height profile from a given dimer covering of the 
hexagonal lattice. Although the dimer covering in turn is 
uniquely determined by the corresponding spin state, we 
have not yet established a direct relation between the spin 
variables (Ji and the displacement field u. We will first 
express the spin variable in terms of the height profile. 
For a given spin state, we chose one up spin and consider 
its site as the point of origin with height h = 0. The con- 
struction of the dimer covering from the spin state means 
that the height on all other sites is then fixed since paral- 
lel (antiparallel) spins imply a height change by +2 (—1) 



from site to site if the down pointing triangles are tra- 
versed clockwise. If a single spin is flipped, the height 
on that site changes by ±3. From this is it clear that 
the height profile can locally be changed by 6 without 
affecting the spin state. However, a spin flip with such a 
change in h leads out of the ground state manifold. Al- 
lowing for this excitations, the height proflle is no longer 
single- valued. This induces topological defects which are 
vortex-anti-vortex pairs in form of two triangles, each 
having three parallel spins, see Fig. [3 Along a closed 
curve around a single (anti-)vortex the height changes 
by 6. One easily proves that the above constraints are 
met by the spin-height relation 



cos q • r i 



±1, 



(25) 



with q = (47r/(3a),0) and the triangular lattice sites 
r;. The relation to the string displacement field is then 
given by the coarse grained relation h = Siu — uq) / i with 
£ — 3a/ 2 in the isotropic case considered here. Using 
these relations, the spin-spin correlation function can be 
written as 



An , 
3a 



{Xi Xj) 



l([u(r.)-u(r,)]=^) 



(26) 



where we have used the fact that u(r) is Gaussian and 
that there are no topological defects at T — 0. From 
Eq. (|24() and ^ciiC44 = tt/P it immediately follows that 



An 
3a 



(27) 



with 77 = 1/2 which is in agreement with the exact 
result''^. The anisotropy factor depends here only on 
the a;-coordinates of the spins since the x-axis had been 
chosen to coincide with one of the triangular lattice di- 
rections. 

Finally, we consider the height correlations for gen- 
eral anisotropic dimer weights. These correlations have 
been studied in the context of the triangular SOS model 
some time agc^^-'^^. In those works, the correlations were 
obtained from four-spin correlations of the TIAF in the 
isotropic case, and in terms of the Pfaffian method for 
the anisotropic situation. In contrast, here we present a 
simple derivation of the correlations based on our effec- 
tive elastic description of the string lattice. The mean 
squared height difference between two distant positions 
increases logarithmically as the string displacement does. 
The coefficient K measures the large-scale stiffness of the 
height profile. 



((Mr)-M0))^) = ^log(9, 



(28) 



By making use of the relation between the height profile 
and the string displacement field, one easily gets 



K 



VC11C44 



9 ' 



(29) 
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(b) 

FIG. 5: (color online) (a) Vortex-anti-vortex pair connected 
by a string loop, (b) Loop connected to a string as to generate 
a non-directed string. This configuration is obtained from the 
one in (a) by flipping a single spin (encircled) which costs no 
energy. 



which is independent of the dimer weights zi and Z2- 
This remarkable universality was found alread y in the 
exact solution of the dimer model in Refs. 122133. but the 
physical reason for that remained unclear. Our interpre- 
tation of the stiffness K as the geometric mean of the 
two elastic constants of the string lattice can explain the 
universality. By changing the dimer weights, one tunes 
both the string density and the string tension g. Naively, 
one would expect that an increasing g would render the 
height profile stiffer. However, there is also a reduction of 
the entropic repulsion between the strings which accom- 
panies the reduction of string fluctuations. Interestingly 
but effects act together as to generate universality since 
C44 ~ ff, Cii 1/g and the density dependent relation 
between h and u. 



D. Finite temperature 

For the TIAF at finite temperatures, the mapping from 
spin states to string lattices can still be applied. The 
strings remain non-crossing, since each triangle can have 
only one or three frustrated bonds. However, the tri- 



angles with three frustrated bonds are excitations that 
generate topological defects, see Fig. Ela). Two defects 
form a pair that spans a string loop which is confined be- 
tween the strings, i.e., strings cannot cross loops. How- 
ever, loops can attach to strings resulting in non-directed 
strings, see Fig.|SJb). The fugacity for the lowest energy 
single-spin excitations is e~^'^l^ . The quasi- long-range- 
order of the T = phase is destroyed if the defect pairs 
can unbind. This will happen at sufficiently weak string 
lattice stiffness. Actually, the critical stiffness for that 
is s/c\\C\l = l-Kji^ and hence by a factor of two larger 
than the actual stiffness corresponding to the TIAF at 
T = 0. Hence topological defects are always a relevant 
perturbation but they can occur only for T > where the 
fugacity is finite, leading to exponentially decaying spin 
correlations. Thus the system does not show a Kosterlitz- 
Thouless transition by tuning the temperature. 

IV. STACKED (2+l)D ISING MODEL 

After we have introduced the description of the frus- 
trated 2D TIAF in terms of fluctuating strings, we shall 
apply this mapping now to study the (2 -I- l)-dimensional 
stacked TIAF which consists of ferromagnetically cou- 
pled TIAF layers, see Eq.lQ). As in the 2D system, the 
in-plane frustration is expected to have strong influence 
on the underlying physics due to the Ising symmetry. 
Actually, the behavior is independent of the interaction 
along the stacking direction, whether it is ferro- or an- 
tiferromagnetic. Experimentally, the stacked model is 
a reasonable description of triangular cobalt antiferro- 
magnets of the type AC0X3 where A is an alkali metal 
and X a halogen atomS. For these magnets, a strong 
crystal field splitting leads to an effective spin 1/2 state 
with the moment oriented along the stacking direction. 
In these compounds the in-plane exchange coupling Jy 
is small compared to the inter-plane coupling J±. The 
increased dimensionality of the stacked system reduces 
the degeneracy of the ground state so that the entropy 
per spin vanishes at T = and ordering at even finite- 
temperature might be possible. This is also suggested by 
the following argument Each ground state of the 2D 
TIAF yields a ground state of the stacked system if all 
spins are aligned along the stacking direction, i.e., the 
configuration in each layer is identical. Hence, at T = 
the spin correlation function cannot decay in-plane faster 
than that of the 2D system. However, one might wonder 
if the system can gain entropy by introducing domain 
walls parallel to the layers which would destroy the order 
along the stacking direction. The average energy cost for 
such a wall is at least 

l-N'/' 27rr 

- 2Jj_ / -^dr - iV^/^ (30) 
Jo 

where N is the total number of spins of the 3D sys- 
tem and we used the power law of Eq. (|27|) . This has 
to be compared to the entropy gain which growth only 
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~ (1/3) logiV since the wall can be placed at any of the 
N^^^ layers. This naive argument suggests the absence of 
domain walls at T = 0, and order along the stacking di- 
rection. Moreover, as we will see more clearly below, the 
huge in-plane degeneracy can actually induce order even 
in-plane, a phenomenon known as order from disordeA 
The stacked TIAF has been studied first by 
Blankschtein et afiS. Their results based on the LGW ap- 
proach and Monte Carlo simulations suggested the exis- 
tence of two different ordered phases and a XY-like tran- 
sition into the paramagnetic phase. The LGW Hamil- 
tonian was constructed for the large scale fluctuations 
about the two minimal energy modes with wave vectors 
Q± = (±47r/3,0), leading to a 3D XY model with a 6- 
fold symmetry breaking term, 

i?LGw(^) = yd3r{i|VVoP + ^|^o|' + U4|Vo|' 

+ U6 1 Vol' +1^6 1^0 1' COS (60)}, (31) 

for the complex order parameter tp{r) = V'oe"^. Since 
for p ^ 3.4 the symmetry breaking term is irrelevant in 
3D |33, the LGW theory predicts for the transition to 
the paramagnetic phase XY universalitj*iSi^. The sign 
of vq is not fixed in this approach, and hence two or- 
dered phases with a relevant symmetry breaking term 
are possible in principle (see below for the two types of 
ordering). In fact, the two corresponding phases were ob- 
served with increa sing temperature in the Monte Carlo 
simulations in Ref. Il6l However, more recent simulations 
indicated the existence of only one ordered phase which 
corresponds to the one found in Ref. at higher temper- 
atures onljii^i2Si2S42L2&. This conclusion is also supported 
by a hard-spin mean field theorji^. There is also a con- 
troversy about the nature of the transition to the para- 
magnetic state. While the simulations of Heinonen et al. 
indicate tricritical behavior—, a histogram Monte Carlo 
analysis supports XY-like behaviopii-. These conflicting 
results have to be viewed against the background of Cop- 
persmith's worlii^. She argued that the above mentioned 
LGW approach is not reliable since it ignores the geomet- 
ric frustration and hence does not yield the correct low 
temperature state. 

In the following, by making use of the string mapping 
established in the preceding section, we provide a mi- 
croscopic derivation of the LGW action for the stacked 
TIAF. This will allow us to discuss the nature of the or- 
dered phase and the transition to the paramagnetic state. 
In each layer, we relate the spin variables aik to a height 
profile hik as in the preceding section so that the relation 
of Eq. (|25|l holds for every layer. Then the intra-layer 
coupling can be written as 



(32) 



with a shift rjij ~ +1 for the bond directions (a, 0) and 
77,;j — —1 for the directions (a/2, ±\/3a/2). Next, we 
have to define how the height should change along the 



stacking direction. Since hik changes by ±3 for a spin 
flip, we set along the column passing through the origin 
of each layer h^k — (3) if a^k = +1 (—1), where k 
numbers the layers. According to this rule, the inter- 
layer coupling reads 



CifcCTifc+i = - COS 



{hik — hik+i) 



(33) 



Now, the height profile can again be interpreted as 
the displacement field of strings, which now form a 3- 
dimensional lattice. Hence, it is expected to more stable 
against fluctuations, and spin order will occur. How- 
ever, the periodic couplings of the displacement field al- 
low for topological defects which eventually drive the sys- 
tem into the paramagnetic state. Even in the presence of 
defects, the picture of non-crossing strings remains valid 
since a triangle can have either 1 or 3 frustrated bonds. 
When introducing continuous fields, the lock-in potential 
of Eq. (|23|l has to be included as to reflect the original dis- 
creteness of the height which reads hik = 3(uifc — uq) /£ in 
terms of the string displacement Uik- Hence, we obtain 
the reduced 3D string Hamiltonian 

i7s = -^11^ cos j{uik -Ujk + ri,ja/2)^ 



K±_'^cos^J{uik-u^k+l) + i^y^cos 



Git 



(34) 



with couplings -fCy = J\\/T^ K± = J±IT and mean string 
separation I — 3a/2. Here v > Q due to the entropy ar- 
gument at the end of Section IIII Bl The in-plane shift 
rjij reflects frustration since X^^y — 3, where the sum 
runs over the bonds of a down pointing triangle. Since 
in the discrete version the Uik can vary over the bonds 
only by +a or —a/2, the energy is minimized for a non- 
uniform change of Uik along the triangles. This is dis- 
tinct from a priori continuous field of an antiferromag- 
netic XY model on a triangular lattice whose oriented 
uniform change defines a helicity, giving rise to an ad- 
ditional Z2 symmetry-8.. Hence, the latter XY model 
exhibits a Ising transition because of the Z2 symmetry 
breaking of the ground states^. However, in the present 
case, the non-crossing restriction for the strings, originat- 
ing from the geometrical frustration, prohibits this extra 
symmetry breaking. 

Thus the stacked Ising model of Eq.(^ maps to a 
stack of planar lattices of non-crossing strings which is 
described by a (2-|-l)D frustrated XY model with a 6- 
fold clock term. Interestingly, this provides a microscopic 
derivation of the LGW theorjii^ if the string displacement 
u is identified with the phase <j) of the order parameter via 
(j) = -Ku/ L However, our effective model differs from the 
LGW theory in two important points: (i) the in-plane 
XY coupling is frustrated and (ii) there is a topological 
constraint on (j) since u is restricted by the non-crossing 
condition. The latter point will lead to an entropicly 
increased phase stiffness on large length scales. 
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(a) 




(b) 

FIG. 6: (color online) Different type of topological defects, (a) 
vortex-anti-vortex pairs in the layers form vortex loops per- 
pendicular to the planes, (b) vortex loop parallel to the planes 
which borders the mismatch regions in which the strings in 
two adjacent planes are shifted by 21. 



With the model of Eq. H34|l at hand, we can analyze the 
ordering mechanism and the transition to the paramag- 
netic phase. The latter transition is driven by topological 
defects which are generated by thermal fluctuations. It is 
well known that the XY coupling allows for defect loops 
across which changes by 27r. These defects are super- 
positions of two kind of loops [cf. Fig.|B]: fully frustrated 
triangles with all spins aligned form in-plane vortex-anti- 
vortex pairs, which can be viewed as a string forming an 
in-plane loop. Together with other pairs in the neighbor- 
ing layers they form defect loops oriented perpendicular 
to the planes, cf. Fig. |HIa). Another type of defect loops 
is oriented parallel to the layers. They arise as the bound- 
aries of 2D areas across which strings in adjacent planes 
are shifted by 2£, cf. Fig.|Hfb). For the universality class 
of the transition to the ordered state it is important to 
take into account the stacked nature of the (2 -I- 1)D XY 
model and the 6-fold clock term. In turns out that the 
stacking is irrelevant since the layers cannot decouple in- 
dependently from the vortex unbinding transition within 
the layers. This is due to the fact that unbinding of de- 
fects oriented parallel to the layers can occur only at a 
critical value for the effective coupling which is by 
factor of 1/8 below that for the in-plane dissociation of 
vortices^. In addition, a 6-fold clock term is known to be 
irrelevant at the XY critical point under renormalization 
if p > 3.4 [s^. Hence we conclude that the transition 
from the paramagnetic phase to an ordered state must 
be in the 3D XY universality class. 



In the following, let us discuss the property of the or- 
dered phase. In the latter phase, the XY couplings of 
Eq. (|34|) can be expanded in Uiu, and in the continuum 
limit each layer is described by Eq. 121|l with an addi- 
tional harmonic inter-layer coupling. Moreover, the 6- 
fold clock is relevant, and the strings are locked into a 
flat state with spin order whose nature depends on the 
form of the locking potential. In order to compare our 
approach with the LGW theorjii^, we discuss potential 
ordered states. As we have seen in Section UlI Bl there 
are two classes of flat strings. They differ by their value 
of the shift uq in the lock-in potential Vl (w) . After going 
over from the discrete field Uik to a coarse-grained de- 
scription with a continuous field u{y), the coarse-grained 
spin variable can be written as 

cr(Rjfe) = cos • Rjfe + ^(M(Rjfe) - uq)^ , (35) 

where Rjfc are the 3D lattice sites and Q — 
(47r/(3a), 0, 0). For flat strings, the coarse-grained fleld 
w(r) is zero on average, and the global shift wq determines 
the sublattice magnetizations. For the first state, cf. left 
part of Fig. ^ uq — —a/ A, and there is an additional shift 
of 7r/6 in Ea. H35() . leading to the sublattice magnetiza- 
tions (cr) — , — ^ , 0) . For the second state, cf . right 
part of Fig. H Mo = 0, one has the sublattice magne- 
tizations (cr) = (1,— 5,— ^) since QR^j, = 0, 47r/3, or 
— 47r/3 on the three sublattices. Hence all sublattices are 
at least partially ordered. Since our entropy argument 
lead to uq = —a/4, the first state with one fully disor- 
dered sublattice should set up the ordered phase at all 
temperatures which appears to be consistent with more 
recent simulations. 

How can this ordering be reconciled with the fact that 
each ground state of the highly degenerate and hence dis- 
ordered 2D TIAF is also a ground state of the stacked 
model when all spins are aligned along the stacking di- 
rection? This is the point where the mechanism of order 
from disorder comes into play. As we have seen at the 
beginning of this section, the formation of domain walls 
parallel to the layers is energetically not favorable. How- 
ever, the system can gain entropy by flipping single (so 
called flippable) spins in a single layer which does not cost 
energy within the layer. Due to the huge ground state de- 
generacy there are many conflgurations which just differ 
in their orientation of flippable spins. Hence, a subset of 
all the 2D ground states is selected which allows for those 
single spin flips. This explains the existence of the found 
ordered state since it is composed only of an appropriate 
subset of all 2D ground states which on average has one 
sublattice disordered. The spins on this disordered sub- 
lattice form chains along the stacking directions which 
are decoupled from each other at low temperatures, and 
hence should behave as individual ID Ising spin chains. 
At sufficiently low temperatures, the stacked system is 
thus expected to show excitations which are ID Ising-like. 
This conclusion is consistent with multi-spin Monte Carlo 
simulations of the specific heat at low temperatures^^. 
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V. QUANTUM FRUSTRATED MODEL 

Having established the relation of the (2 + 1)- 
dimensional TIAF to a lattice of fluctuating strings with 
topological defects, we now will apply this relation to 
the 2-dimensional TIAF in a transverse magnetic field. 
The Hamiltonian of this model is given in Eq. ||2J). The 
transverse field introduces simple quantum dynamics to 
the highly degenerate critical state of the classical TIAF. 
Due to the general relation between 2D quantum spin 
systems and 3D classical Ising spin models^^, one can 
expect from the results for the stacked TIAF quantum 
order arising from the interplay of quantum fluctuations 
and geometric frustration. The existence of a quantum 
ordered state was suggested in recent works on the basis 
of a LGW theory?-, a kink modeliS, and simulationgi. 

The exact mapping from the 2D quantum TIAF to 
a classical stacked Ising system is provided by Suzuki's 
theorem which states that the partition function of the 
quantum system can be written asii 



Z = 



Trexp<'-^5]afa| + ^;^af 









lim 


- sinh 




n — >oo 
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\nT ) _ 



ilogCOth j ^<y^,k(y^,k^'^, 



(36) 



where (Tik = ±1 are classical Ising variables which are 
defined on a stacked system of n antiferromagnetic tri- 
angular lattices, and the second index k = 1, . . . , n refers 
to the layer number. Hence, in the relevant limit of large 
n the quantum system is described by the Hamiltonian 
of the stacked TIAF of Eq. 0J with the couphngs J|| 
and Jl replaced by a infinitesimal small in-plane cou- 
pling i^ii = J/{nT) and a strong ferromagnetic inter- 
layer coupling K± — i log(rtr/r) which is controlled by 
the strength of the quantum fluctuations, i.e., the trans- 
verse field r. In analogy with the preceding section, the 
quantum system is hence described by a (2 + 1)D string 
lattice with the Hamiltonian of Eq. H34I) . This suggests 
that, as for the classical 3D system, the transverse field 
Ising model must have an ordered phase at small trans- 
verse field and a quantum phase transition in the 3D XY 
universality class to a disordered state. This scenario was 
previously predicted from by the LGW approach. How- 
ever, in contrast to the latter approach, our mapping of 
the quantum frustrated spin system to strings yields even 
the coupling constants of Eq. H34II which is of utmost im- 
portance for the construction of the phase diagram. 



A. Quantum critical point 

First, let us focus on the quantum phase transition to 
the paramagnetic phase. The purpose of the following 



analysis is to estimate the location of the critical point 
and to show that the Trotter limit ti — > oo can be carried 
out explicitly in our approach. It is useful to separate 
the partition function associated with the XY model of 
Eq. H34|l into a spin wave part and a vortex part. This 
can be done by making use of the fact that the planar 
model can be mapped to a simpler model proposed by 
Villain^^. As was shown by Kleinert^^ and in Ref. 45, in 
the partition function one can make the substitution 



^Kcos{4>) 



oo 
m— — oo 



-K{(i)-27zmY /2 



(37) 



where the right hand side is known as Villain 
couplingMii^ with a new iiT-dependent coupling K . Both 
expressions become identical in the two limits of ^ 
and if — > oo which are fortunately precisely the two cases 
arising from the Trotter limit n — > oo. Then the coupling 
constant K and the coefficient c are given by 



K 



K = K, 
1 



c = for .R" oo 



(38) 



21og(2/X)- 



y/2T:K for K ^0 (39) 



Hence, after taking the Trotter limit, the coupling con- 
stants of the equivalent Villain model both scale logarith- 
mically with nT, 



Kn 



1 



2 log(2nr/ J) 



(40) 



Only for T = 0, one can set T ^ J/n so that for n ^ oo 
the coupling K± remains finite, and the system shows 3D 
behavior. For any finite T, however, K± must diverge 
with n and the system shows a 3D to 2D crossover with 
increasing length scales. 

Knowing the coupling constants of the Villain model 
with decoupled spin wave and vortex parts for large n, we 
can go ahead and apply dimensional crossover scaling in 
order to make quantitative predictions about the location 
and universality of the quantum critical point. We start 
by eliminating n from Eq. (|40|l . From this we obtain the 
temperature independent relation 



4^ \ 



(41) 



between i^Ty and K±^ which makes the parameter space 
one dimensional. The geometric mean of the coupling 
constants controls the stiffness of the XY system. It can 
be expanded in the relevant limit of small Kn as 



1 - log(2r/ J) K« + 0{Kl 



(42) 



where we used the relation of Eq. H41() . The latter rela- 
tion will be compared to a crossover scaling result which 
provides a quantitative description of the increase of the 
transition temperature of a layered (2-t-l)D XY model 
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under the increase of the number of layers. Following 
closely the analysis in Refs. i4(ai47>. one obtains a relation 
between the 3D critical value for the in-plane cou- 
pling and the corresponding if^ for a system of n layers, 



1 K"" 
n KS. 



7 



K± 



1/2 



1 - 



(X 



(43) 



with the critical exponent « 2/3 of the 3D XY model 
and a numerical constant 7. For n = 1 this relation 
yields at the 3D critical point with = the fol- 
lowing expression for the geometric mean of the coupling 
constants, 




(44) 



where we have expanded for <C Kf. For consistency, 
the latter expression should be of the same form as the 
expression in Eq. (|42|l at the quantum critical point with 
r = Fc corresponding to ify — K^. Indeed, both ex- 
pressions are identical if one sets 



7=l/(2if^), 2r,/J = e''/^^. 



(45) 



Hence, the location of the quantum critical point is de- 
termined by the critical value for the coupling of the 
2D XY model on a triangular lattice. The standard 
Kosterlitz-Thouless argument for the vortex unbinding 
transition yields'*^ in that case = 2/7rx2/V3, leading 
to Fc/J = 1.24. However, this scaling approach neglects 
renormalization effects due to the 6-fold clock term and 
the non-crossing of strings which should provide a net in- 
crease of Fc. This is consistent with recent Monte Carlo 
studies which suggest Fc/J ~ 1.65 0- 



B. Phase diagram 

Having established the existence of an ordered phase 
for F < Fc at T = 0, it is important to study the sta- 
bility of this phase against thermal fluctuations. For fi- 
nite temperatures, 2D XY physics should dominate at 
large length scales, and two Kosterlitz-Thouless transi- 
tions separating a critical phase at intermediate temper- 
atures from the ordered and the paramagnetic state, re- 
spectively, are expected^. In the following we will con- 
struct a quantitative phase diagram for the transverse 
field TIAF in the F-T-plane. The stabihty range of the 
phase with bound defects can be estimated from the 
crossover scaling formula of Eq. H43(l . Setting in the lat- 
ter expression K^^ to its critical value for a n-layer 
system, and using the scaling formula also at n = 1 in 
order to relate to Kf and the bare coupling ratio 
K±/K\\, we can derive an expression for the Kosterlitz- 
Thouless transition temperature Tc. 2(F) of the critical 



FIG. 7: (color online) Phase diagram with the phase bound- 
aries given by Eq. 1461 with b — 0.98 and Monte Carlo results 
of Fig. 1 in Ref. 7. A critical phase (C) is separated by a 
boundary at Tc,2 from the paramagnetic phase (PM) and at 
Tc.i — 4/9 Tc. 2 from an ordered phase (O). 



phase (C) which has only bound defects [see Appendix 
^for details], 



J Tc ^ \ r 



(46) 



where & is a numerical constant which is fixed by the 
(unknown) renormalization of K± and remains finite for 
n — i- 00. At large but finite n, for consistency, the 
renormalized effective coupling for the n-layer system 
must behave at as ~ {n/ric)^/^ with a characteris- 
tic number Uc = bKf{J/T){T/Tc) which is a measure for 
the strength of quantum fluctuations and characterizes 
the effective system size along the Trotter ("imaginary 
time") axis. The spin correlation function in the critical 
phase (C) decays according to Eq. H27|l with the expo- 
nent ry varying continuously between 77 = 1/4 at Tc_2 and 
7y = 1/9 at a lower critical temperature Td which marks 
the transition to the ordered phase (O)--. The correla- 
tion function exponent behaves discontinuously at F = 
where it is ?/ = 1/2 and at the quantum critical point 
where the 3D XY resuh 77 « 0.040 holdaiS.. 

At the lower critical temperature Tc.i there is a second 
Kosterlitz-Thouless transition to an ordered state where 
the clock-term in Eq. H34|) is relevant and locks the strings 
to the lattice. For a single layer, this transition occurs 
at the critical couphng = 9/(27r) x 2/^3, see Ref.llsl 
The boundary of the ordered phase can be obtained anal- 
ogously to that at Tc,2, yielding 



TcA 



9 



(47) 



Close to the quantum critical point both Kosterlitz- 
Thouless transition temperatures vanish ~ (Fc — F)'' 
as expected from scaling. Fig. \7\ compares our analyt- 
ical results to recent Monte Carlo data for the phase 
boundaries, showing very good agreement across the en- 
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tire range of F if we set b = 0.98 in Ea. (|46|l . Accord- 
ing to our analysis of the ordered phase for the stacked 
system in Section IIVI the quantum ordered phase (O) 
is characterized by the finite sublattice magnetizations 
(\/3/2, — \/3/2, 0). This type of order is consistent with 
recent simulations^. 



VI. CONCLUSION AND EXTENSIONS 

In this article, we have studied for a 2D Ising system 
the combined effect of classic geometric frustration and 
an ordering mechanism which either results from stack- 
ing many 2D system or by allowing for quantum dynam- 
ics from a transverse field. We have presented an ex- 
act relation between spin variables and the displacement 
field of a string lattice which we used to derive an effec- 
tive Hamiltonian which was obtained previously only on 
the basis of a LGW approach. This allowed us to de- 
termine the nature of the ordered phase, the universal- 
ity of the quantum phase transition to the paramagnetic 
phase, and the quantitative phase diagram at finite tem- 
peratures. To our knowledge, this is the first approach 
for the analyzed models which explicitly performs the 
Trotter limit rt — > cxd to obtain quantitative results for 
the transition temperatures. We have compared the so 
obtained phase diagram to predictions of recent simu- 
lations and found satisfying agreement. We related the 
spin model also to its dimer model representation and 
the resulting height profile description. Implementing a 
simple random walk argument for the strings, we could 
derive the exact free energy in the continuum limit which 
was known before only from a more complicated Pfaf- 
fian method. Our approach could also explain the phys- 
ical reason for the previously found universality of the 
large-scale stiffness of the height profile in the presence 
of anisotropic bond weights. 

For future extensions of our approach it is interesting 
to note that the string analogy applies independently of 
the form of the spin couplings. Interestingly, the map- 
ping to strings should also be useful in the presence 
of quenched disorder when the system usually exhibits 
spin glass behavior. The behavior of a 2D string lat- 
tice with pinning site^ suggests that in geometrically 
frustrated syste ms the sp in glass correlation function^ 
gsG (ri — Tj) = < aiCFj >2 decays at least according to a 
power law, which implies that the spin glass phase is de- 
stroyed by the huge degeneracy from the geometric frus- 
tration. However, quenched disorder will induce topolog- 
ical defects, and a full analysis is more complicated. The 
situation is under better control for geometrically frus- 
trated Ising systems with random dilution. The dilution 
sites act as pinning centers for the strings and topological 
defects do not arise from this type of disorder. 

In the quantum frustrated system, though it is very 
hard to get a clear picture of how the transverse field and 
quenched disorder act together, one can try to make use 
of a method analogous with what Moessner et alA used 



in discussing the existence of an ordered phase in clean 
quantum frustrated Ising systems. Moessner et al. used 
the divergence of the integrated squared two-spin corre- 
lation function to argue about a non-analytic free energy. 
Instead, with quenched disorder, the spin glass correla- 
tion function had to be integrated to detect a singularity 
in the disorder averaged free energy. The divergence of 
the latter integral would then signal the possibility of an 
(spin glass) ordered phase, induced by quantum fluctua- 
tions. 

Finally, we discuss potential implications for the quan- 
tum phase transition in 2D Ising spin glass models in 
a transverse field^. From the analysis in Section [V] 
we know that the 2D quantum frustrated Ising model 
can be mapped to a system of classical stacked string 
lattices which is described by a 3D XY Hamiltonian. 
With quenched randomness, the quantum spin system 
should be mapped to a 3D XY model with columnar, 
i.e., correlated disorder. If we assume that the lattice 
type (and hence the geometric frustration) is irrelevant 
for the universality class of the transition, we can ar- 
gue that the quantum phase transition of 2D transverse 
field Ising spin glasses should belong to the universality 
class of 3D XY model with columnar disorder. Interest- 
ingly, this speculation seems to be consistent with recent 
numerical results for both a transverse field Ising spin 
glassi^ and a 3D XY model with columnar disorder—. 
The two corresponding exponents, the dynamical expo- 
nent z = 1.5 ± 0.05 for the spin glass and the anisotropy 
exponent ^ = 1-3 ± 0.1 for the XY model, are indeed in 
some agreement. 
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APPENDIX A: PHASE BOUNDARY 

In this appendix we derive the boundary of the critical 
phase (C). We start from the scaling result of Eq. H43(l 
with n — 1. This allows us to compute a relation between 
the critical couplings Kf and . In the limit ^ Kf 
we find 

^ -i^JtT T^' (Al) 

where v is the critical exponent of the 3D XY model, and 
7 = 1/(2/^5^). Using the latter result again in the scaling 
formula of Eq. (|43|l for arbitrary n and with X|| = 
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set to its critical value, we obtain 



1 




(A2) 



Multiplying with T/J and using the relation 2Tc = 
j^vjKi i-j-^g quantum critical point, we get the final 
result 



which yields as a function of the inter-layer coupling 
at large but finite n. We continue by introducing the 
bare coupling constant Co — e^^^ so that n — C^oT/T. 
The renormalized effective coupling constant will be de- 
noted by C = e^^-L . This notation allows us to take into 
account the fact that the renormalization of the couplings 
K\\ and K±^ is dependent due to their relation by Eq. (|41|l . 
Then the effective couplings can be written as 



1 



21og(2Cr/J)' 



if*? 



^logC. 



(A3) 



Now we use = K^^ and = Kf in Eq. and 
expand the resulting equation for large C which gives 



1 = 



Co 



r 



{2\ogQf/^ KIT \K1 ^°^(^J 



(A4) 



J 



Coe''/^^ 



2i^f(21ogC)5/3 r. 



log^ 



(A5) 



which corresponds to Eq. (|46|l if we denote by h the coef- 
ficient in the latter expression. For consistency, it follows 
from = i log C and the definition of h that the ef- 
fective inter-layer coupling must diverge with n — > c» 
according to K'^^ = [n/ricf'/^ with a characteristic 
Tie = bKfJT/{TTc) which itself diverges for T ^ 0. 
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